Current driven spin- wave instability triggered by the anomalous Hall effect. 
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We studied the effect of strong electric current on spin waves interacting relativistically with the 
current. The spin- wave spectrum is calculated at arbitrary direction of the wave vector. It is shown 
that the alternating Hall current generated by the alternating magnetic moment of the spin waves, 
reduces the spin-wave damping. At strong enough unpolarized dc current the damping changes sign, 
and the spin-wave amplitude starts to increase exponentially fast with time. The critical current 
for the spin-wave instability is determined mainly by the anomalous Hall effect, and can be much 
smaller than that for the spin-torque mechanism of instability. 
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Current induced switching of the magnetization or spin 
wave excitation in magnetically inhomogeneous systems, 
e.g. multilayers, has received considerable attention dur- 
ing the last decade. [H When a spin-polarized current 
passes through a ferromagnetic layer, the electrons trans- 
fer their spin angular momentum to the localized spins 
of the ferromagnet resulting in a spin torque acting on 
the magnetization. 0, Q The theoretical predictions were 
confirmed by several groups, the experiments being per- 
formed mainly with trilayers of structure ferromagnet- 
normal metal-ferromagnet, in which the layer magneti- 
zation may be noncollinear, see Ref. rH and references 
therein. 

It has been argued recently that polarized current can 
affect the magnetic properties also of an homogeneous 
bulk ferromagnetic metal. HUSH Adding the spin 
torque linear in the spin current, to the Landau-Lifshitz 
equation of motion, one gets a modified spin-wave (SW) 
spectrum, which shows a current driven instability. In a 
half-metal, when the density of the minority carriers is 
zero, the spin current is equal to the electric current. In 
this case the uniform ferromagnetic state becomes unsta- 
ble at a critical current given by the relation 0, [B| 



Vd = LUk, 



(1) 



where k is the spin- wave wave vector, is the spin- wave 
dispersion, and Vd is the electron drift velocity, which is 
proportional to the electric current. This "Doppler-shift" 
critical current is of order of 10 9 A/cm 2 ,[5( but in general, 
as has been shown by Tserkovnyak et aL,[7[ the critical 
current can be strongly enhanced. 

High enough current densities can excite SW exci- 
tations in ferromagnetic layers even when the current 
is unpolarized, if the source and drain contacts are 
nonsymmetric.jH SW excitation by an unpolarized cur- 
rent injected into a single ferromagnetic film from a point 
contact was observed by Ji et al. [9( and considered the- 
oretically in Refs. 

The original models by Slonczevski Q and Berger Q 
and all subsequent considerations of current induced SW 



excitation rely on the exchange model of interaction be- 
tween the itinerant electrons and the localized spins. In 
this paper we concentrate on the electromagnetic (rela- 
tivistic) interaction of the electron current with the field 
of the SW, which does not preserve the total spin. We 
show that in ferromagnetic conductors with large anoma- 
lous Hall effect this interaction can lead to current in- 
duced SW instability at critical unpolarized current of 
the same order or even smaller than that in the exchange 
interaction models. Unlike the exchange coupling of itin- 
erant electrons with the SW, which is effective only in 
the vicinity of the interface between normal and ferro- 
magnetic layers, |3J the above relativistic interaction acts 
also in the bulk of the ferromagnet. Therefore, the cur- 
rent induced SW instability cased by the Hall effect is 
not restricted to a layered structure of the ferromagnet. 

The physical mechanism for SW generation via the rel- 
ativistic interaction is as follows. Suppose for simplicity 
that an electric current jo is driven parallel to the mag- 
netization Mo. Consider a spin wave propagating along 
Mq. The oscillating magnetic moment and magnetic field 
of the spin wave lie in the plane perpendicular to jo and 
Mo. This gives rise to an alternating Hall current per- 
pendicular to Mo, which in his turn creates a magnetic 
field amplifying the field of the wave. If the electric field 
is strong enough, the amplification will exceed the damp- 
ing due to eddy currents. If there are no other sources of 
damping, the spin-wave system becomes unstable at such 
field. We show in this paper that the instability is not 
restricted to the above simple geometry. It takes place 
at any mutual orientation of jo and Mo . 

In ferromagnetic metals the main contribution to the 
Hall current comes from the anomalous Hall effect caused 
by the spin-orbit coupling in the metal. Since the anoma- 
lous Hall constant is by orders of magnitude larger than 
the normal one, the critical current for SW amplifying 
can be relatively small. 

The full set of equations, which describes SW in a con- 
ducting media interacting with an electric current con- 
sists of the Landau-Lifsitz equation and of the Maxwell's 



2 



equations. The Landau-Lifshitz equation for the preces- 
sion of the magnetic moment is|ll| 

9m lT 1 1 

_= 7 MoXH e// + -H e// - — 



2M x(M xH e// ) 



(2) 

Here m is the (small) transverse alternative part of the 
magnetization, m <C Mo, 7 is the gyromagnetic ratio, r\ 
and T2 are phenomenological SW relaxation times, and 
H e / f is the effective magnetic held given by 



H 



eff 



D 



m • H 



7M0 M 
K(mo ■ n) 2 m + Kn(n ■ m), 



(3) 



where h is the alternative part of the magnetic held, D is 
the stiffness constant of the SW, Ho is the external mag- 
netic held, which in restricted samples includes also the 
demagnetizing helds, n and mo are unit vectors directed 
along the anisotropy axis and the magnetization respec- 
tively, and K is the dimensionless anisotropy constant. 
The Maxwell equations are: 

_ , 47T. 

V x h - — j, 

c 

19b 

Vxe = 

c at 

Vb = 0. (4) 

Here the alternating magnetic induction is b = h + 47rm, 
e is the alternating electric held, c is the light velocity, 
and j is the alternating electric current given by [l2| 



j = ere + a(R B 3o x b + R M 3o X m), 



(5) 



where Rb and Rm are the ordinary and anomalous Hall 
coefficients respectively, a is the conductivity, and jo 
is the dc part of the electric current density. We ne- 
glected in Eqs. Q the displacement current, since for 
conductors considered here the inequality lo -C cr always 
holds. We skipped in Eq. |(SJ) a small term of order 
<j(RbBq + RmMo) -c 1. The SW frequences considered 
in the paper are of order of (10 9 — 10 10 )s _1 . i.e. much 
smaller than all typical electron relaxation frequencies in 
a ferromagnetic metal. Therefore, we used below the dc 
values of the transport coefficients. 

The Maxwell equations, with the current from Eq. J5J), 
relate the Fourier transforms of h and m as 
2 

h = -47r[w fc 2 + « — (k • v 6 - w)] _1 [£j (k • m)k 



S 2 
2 

-i — (k • m)v„ 



. 2 



kv)m]. (6) 



Here loq = ^ttjMq, the skin penetration depth, S, at fre- 
quency ujq is given by S = (c 2 jl-Kauj^) 1 ! 2 , and the ef- 
fective velocities v&,v m and v are related to the Hall 
coefficients by 



v fa = R B Qo = v d , v„ 



RmcJc 

47T 



(7) 



It is supposed that S is much larger than the electron 
mean free path. 

The Landau-Lifshitz equation yield another relation 
between h and m. When Ho and Mo are parallel and 
directed along the anisotropy axis, this relation reads: 

(f^k — lo — iaD.\ s )m,+ — 7_Mo/i+(l — ia) 

(fi k + oj + iaQk)m- = •yM h-(l+ ia). (8) 



Here 



1 



1 



1 



7-Mq V 1 "! T 2 



(9) 



m± = m x ± iniy , h± = h x ± ih y , while the axis z is along 
the magnetization. The frequency Qk is given by 



n k = 7 (H Q + H a ) + Dk 2 , 



(10) 



where H a is the anisotropy field: H a = KMq. 

Eqs. ([!]) and © give the dispersion relation for SW 
in an external electric held. We consider in what fol- 
lows wave- vectors k, which are larger than S^ 1 : M > 1. 
One obtains then to leading order in the small parameter 
l/fc 2 £ 2 : 

Wk(jo) = ^k^ik + -^jpr( k x v m) • m 
■(fii k + fi k cos 2 e)(wt(0) -k-v) 
+acj k (0)(2fi k + w sin 2 9) 



5 2 k- 



flkk ■ m [(v m • k)(m • k) - fc 2 v TO • m ]. (11) 



Here 9 is the angle between the wave vector k and the 
magnetization, Oik = Ok+^o sin 2 9, and ^k(O) = o->k(c = 
0) = v^k^ik- In the absence of the current jo, this 
equation gives the usual spectrum of SW decaying due 
to Landau-Lifshitz-Gilbert damping a, and due to the 
eddy currents. The last decay is proportional to l/k 2 S 2 . 

In what follows we consider such wave vectors that the 
contribution to the SW damping from the last term in 
Eq. (fTTj) is equal to zero. This happens specifically, if k 
is along or perpendicular to the magnetization, or at any 
9 provided k is along the current. Eq. (jlip then yields: 

tj f\ / or. 1 ( k x Vm ) ' m " 
ttewkUo) = V"k"ikH ^5 , 



Imcjk(jo) = -(fk + otk) <x»k(0) - 



-k-v (12) 



Here 



fk 

Ok 



Oik + O k cos 2 9 



k 2 8 2 V^Jh^ 

a (2Q k 
2\/OkOik 



ujq sm 



(13) 
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Eq. (fT2|) is our main result. It shows that Imw changes 
sign at a critical velocity, v c , given by the relation 



k-v. = ■ 1 - ^ ) ■ 



(14) 



At higher effective velocities, i.e. at higher currents the 
amplitude of spin waves, with k satisfying Eq IT4")) , in- 
creases exponentially with time. 

When Rm /47r is smaller than Rb, i.e. v — v& = Vd, the 
SW instability condition (fl"4f resembles those, obtained 
in Refs. 0, [1,0, for SW instability in half-metals trig- 
gered by spin-transfer torques, see Eq. ([IJ. However, 
Rm in ferromagnetic conductors is usually by many or- 
ders of magnitude larger than Rb- We have, therefore, 
v « v m 3> Vd, and the critical current for the instability 
considered here is much smaller than the critical current 
given by Eq ((T|). 

Measurements [l3| performed on Fe, Co, Ni and Gd 
films with the thickness of 1 yum, show that in pure met- 
als, with the resistivity p = (10 -4 — 10~ 6 ) Slcm, the 
anomalous Hall conductivity <th is of order 10 3 (f2cm) _1 . 
Thus, for metals with the resistivity p — (10 -4 — 
10~ 5 )fkim, and with M « 10 3 G, one gets R M = 
(10" 8 - 10" 10 ) Qcm/G. The values of R M for Ni films, 
which follow from the data obtained in Ref. [lJ, also fall 
in this region. It follows than from Eq ([7} that the effec- 
tive velocity v m is of order v m = (10 -1 — 10 -3 )jo- Here 
jo is in A/cm 2 , and v m in cm/sec. Note that the typical 
drift velocity is of order 10 -4 jo cm/sec. Q 

The real part of the SW frequency, Eq. (|12p. also 
acquires a term linear in the current, which is solely 
caused by the anomalous Hall effect. The frequency of 
spin waves, with k non-parallel to the magnetization and 
to the current, is modified by the current. The current 
increases or decreases the frequency, depending on the 
direction of k. 

When the wave- vector k and the current are along the 
magnetization, the critical velocity is given by 



v c {k) 



Dk + 



j(H + H a ) 



k 



(i + f«v) 



(15) 



v c (k) is minimal at k — ko given by 
1 



6aS 2 



(2 + ap) + V4 + 28a/? + (a/3) 2 ] , (16) 



where (3 = j(H + H a )5 2 /D. 

When the damping a is small, a <C 1//3, the critical 
velocity coincides with the phase velocity of the SW,[l2| 
and k is given by fc = [j{H + H a )/D} 1 / 2 , v c (k ) being 
equal to 



v c {k ) = 2y/D-y(H + H a ). 



(17) 



With typical values D = (0.1 — 0.05) cm 2 /sec and 
-/(H + H a ) = 2 • 10 9 sec" 1 , one gets v c (k ) « 3 • 10 4 



>° 15 




FIG. 1: Dependence of the minimal critical velocity on a at 
different conductivities a = (1,2,5,10) x 10 4 (ficm) _1 , with 
parameters D — 0.1 cm 2 /sec, Ho + H a = 100 Oe. Both k and 
jo are along the magnetization. 



cm/sec. Thus, the minimal critical current density is 
3c = Anv c {k )/R M c w (3 • 10 5 - 3 • 10 7 )A/cm 2 . This 
value of jc is by several orders of magnitude smaller than 
that obtained in Refs d and @. 

The wave vector fco decreases with increase of a. If a is 
large, a ^> 1//3, it follows from Eq. (Tl7jj) . that ko does not 
depend on D and 7-ffoj and is equal to: k = (2/aJ 2 ) 1 / 2 . 
The critical velocity in this case is 



vM = S 1 (H + H a )V2a 



(18) 



increasing linearly with the external magnetic field. 

With a = 10 5 Om" 1 - cm -1 , and ujq — 10 11 sec -1 , one 
gets 5 = 10 -4 cm. Then, with the above values of D 
and -f(H + H a ) one finds that the inequality a» 1/0 is 
fulfilled, if a is larger than 10~ 2 . The dependence of the 
minimal critical velocity on a at different values of the 
conductivity is shown in Fig. 1. 

If k is parallel to the current and perpendicular to the 
magnetization, Eqs. (fH)) and (fT3"|) yield 



v c {k) = -y^VH^)(l + 



ad 2 k 2 2fl k + 



u> 



Q k +uj 



(19) 



Usually the inequality jHq -C wq holds. Then, at small 
a , a 1//3, the minimal critical velocity is 



z(ko) = y Dujq ~ 10 5 cm/sec, 



(20) 



while the minimal critical current is of order (10 6 — 10 8 ) 
A/cm 2 . As before, the critical velocity increases with 
increase of a. When a is large, and satisfies the inequality 
a ^> 1/(3, the minimal critical velocity is 



v c (k ) = 26 v / 7(#o + H a )uj Q a, 



(21) 
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The critical velocity and the critical current can be con- 
siderably smaller than the above values, if the ferromag- 
net is in a state close to an orientational phase transition 
caused by an external magnetic field. We consider now 
the instability condition in several different arrangements 
of this type. In all cases we suppose that k is along the 
dc current, since we are interested in the minimal critical 
current. 

First, let in a uniaxial ferromagnet the external mag- 
netic field is aligned perpendicular to the easy axis, and 
Hq is larger but close to H a . Then, the magnetization 
points along H . Repeating the previous calculation for 
this orientation of the field, one gets for k and jo along 
the magnetization: 



«•<*>- t( 



1 



where II ill 



= V[Dk 2 + 7 (ffo - H a )]{Dk 2 + jH ) 



At small a 



a < 



D 



7 \H -H a \6 2: 



(22) 



(23) 



(24) 



and at Hq — H a <C H a , the critical velocity is minimum 
at a wave vector, given by 



j 2 H a (H -H a ) 



D{D + a5 2 1 H a ) 



1/4 



and v c (ko) is 



s (fco) = VD-fHa, 



(25) 



(26) 



which is considerably smaller than the critical velocity 
([T7]). At large a: a > D/-f(H - H a )S 2 , the critical 
velocity equals to 



= l6y/2aH a {H -H a ). 



(27) 



Note that if (Hq — H a )/H a is small, the inequalities Eq. 

and kyS 2 <g; 1 are fulfilled whenever a is smaller 
than 1, i.e. the SW damping in this case almost does not 
affect the critical current. 

Suppose now that the ferromagnet is in a metastable 
state, with the field Ho smaller than H a and opposite 
in direction to the magnetization. It follows than from 
Eq. ([3|) that in Eq. fp~5|) and in all subsequent equations 
for v c the field Hq should be replaced by — Hq. Hence, 
when Hq approaches H a the critical wave-vector and v c 
tend to zero, ko is restricted from below by the inequality 
koS 2 ^> 1. This gives for small damping a, which satisfies 
the inequality j(H a - Hq) > DcT 2 , v c > 2D/S ps 

10 3 cm/sec, and j c ^> 10 4 A/cm 2 . This implies that 
a relatively small current of order or larger than (3 — 
5) x 10 4 A/cm 2 can drive the magnetization switching at 



magnetic fields smaller than H a , if the current flows along 
the magnetization. Note that, as in the previous case, the 
inequality Eq. (|24|) is equivalent to the condition q<1. 

Finally, consider thin films, with kd <C 1, d is the film 
width. The SW spectrum in this case at different direc- 
tions of the magnetization and external magnetic field 
was derived in many papers, see e.g. Refs. Ila . We con- 
sider the case, when the external magnetic field is perpen- 
dicular to the film plane and larger than AttMq+H s — H a , 
where H a is the volume anisotropy field, the easy axis 
being in a symmetry direction of the film, and H s is the 
surface anisotropy field. Then, the magnetization is also 
perpendicular to the film, and the spectrum of SW with 
k in the plane is: 



UJ k 



[y(H + 2irM kd) + Dk 2 } [yH + Dk 2 } 



(28) 



where H = Hq + H a - H s - 4ttM . 

As argued above, the damping can be neglected if the 
ferromagnet is in the vicinity of the phase transition. 
Then the critical velocity is equal to the SW phase ve- 



locity, and is minimal at fco = "JlH/D, while H should 
be larger than D/^/S 2 . The minimal critical velocity is 
given by: 

v c (k )= U-yHD + uj d^jHD\ ' . (29) 

The above inequalities yield that v c (ko) is restricted from 
below as: v c > \JTooDdfb w 10 4 cm/sec. Hence, the 
critical current is larger than (10 5 — 10 7 ) A/cm 2 . 

In conclusion, we have calculated the effect of an elec- 
tric current on the SW spectrum in a ferromagnetic 
metal. We have shown that the ordinary and anomalous 
Hall currents lead to the reduction of the SW damping, 
caused by the eddy currents. At sufficiently strong cur- 
rents the damping changes sign, and a SW instability 
develops. The critical current of the instability is deter- 
mined mainly by the anomalous Hall effect, and may be 
much smaller than the critical current for SW excitation 
with spin-transfer torques. 
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